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ABSTRACT 

We review the complications involved in the conversion of stellar luminosities into 
masses and apply a range of mass-to-luminosity relations to our Hubble Space Tele- 
scope observations of the young LMC star clusters NGC 1805 and NGC 1818. 
Both the radial dependence of the mass function (MF) and the dependence of the 
cluster core radii on mass indicate clear mass segregation in both clusters at radii 
r < 20 — 30", for masses in excess of ~ 1.6 — 2.5Mq. This result does not depend on 
the mass range used to fit the slopes or the metallicity assumed. It is clear that the 
cluster MFs, at any radius, are not simple power laws. 

The global and the annular MFs near the core radii appear to be characterised by sim- 
ilar slopes in the mass range (—0.15 < logm/M0 < 0.85), the MFs beyond r > 30" 
have significantly steeper slopes. 

We estimate that while the NGC 1818 cluster core is between ~ 5 and ~ 30 crossing 
times old, the core of NGC 1805 is likely < 3 — 4 crossing times old. However, since 
strong mass segregation is observed out to ^ 6i?corc and ^ 3i?corc in NGC 1805 and 
NGC 1818, respectively, it is most likely that significant primordial mass segregation 
was present in both clusters, particularly in NGC 1805. 

Key words: stars: luminosity function, mass function - galaxies: star clusters - 
Magellanic Clouds - globular clusters: individual: NGC 1805, NGC 1818 



1 PRIMORDIAL VERSUS DYNAMICAL MASS 
SEGREGATION 

The effects of mass segregation in star clusters, with the 
more massive stars being more centrally concentrated than 
the lower-mass stars, clearly complicates the interpretation 
of an observed luminosity function (LF) at a given position 
within a star cluster in terms of its initial mass function 
(IMF). Without reliable corrections for the effects of mass 
segregation, hence for the structure and dynamical evolution 
of the cluster, it is impossible to obtain a realistic global 
cluster LF. 



1.1 Dynamical Evolution in Star Cluster Cores 

Dynamical evolution in dense stellar systems, such as Galac- 
tic globular clusters (GCs) and rich Large Magellanic Cloud 
(LMC) star clusters, drives the systems towards energy 
equipartition, in which the lower-mass stars will attain 
higher velocities and therefore occupy larger orbits. 



Consequently, the high-mass stars will gradually sink 
towards the bottom of the cluster potential, i.e., the cluster 
centre (cf. Spitzer & Hart 1971), with the highest-mass stars 
and those closest to the cluster centre sinking the fastest, 
although this process is not negligible even at the cluster's 
edge (e.g., Chernoff & Weinberg 1990, Hunter et al. 1995). 
This leads to a more centrally concentrated high-mass com- 
ponent compared to the lower-mass stellar population, and 
thus to dynamical mass segregation. 

The time-scale for the onset of significant dynamical 
mass segregation is comparable to the cluster's dynamical 
relaxation time (Spitzer & ShuU 1975, Inagaki & Saslaw 
1985, Bonnell & Davies 1998, Elson et al. 1998). A cluster's 
characteristic time-scale is may be taken to be its half-mass 
(or median) relaxation time, i.e., the relaxation time at the 
mean density for the inner half of the cluster mass for cluster 
stars with stellar velocity dispersions characteristic for the 
cluster as a whole (Spitzer & Hart 1971, Lightman & Shapiro 
1978, Meylan 1987, Malumuth & Heap 1994, Brandl et al. 
1996), and can be written as (Meylan 1987): 
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where -Rh is the half-mass (median) radius (in pc) , Aftot the 
total cluster mass, and (m) the typical mass of a cluster star 
(both masses in Mq). 

Although the half-mass relaxation time characterises 
the dynamical evolution of a cluster as a whole, significant 
differences are expected locally within the cluster. From Eq. 

it follows immediately that the relaxation time-scale will 
be shorter for higher-mass stars (greater (m)) than for their 
lower-mass companions; numerical simulations of realistic 
clusters confirm this picture (e.g., Aarseth & Ifeggie 1998, 
see also Hunter et al. 1995, Kontizas et al. 1998). From this 
argument it follows that dynamical mass segregation will 
also be most rapid where the local relaxation time is short- 
est, i.e., near the cluster centre (cf. Fischer et al. 1998, Hil- 
lenbrand & Hartmann 1998). The relaxation time in the core 
can be written as (Meylan 1987): 

t,,o = (1.55 X 10') , ,, yr, (2) 

(mo) log(0.5 Mtot/(m)) 

where -Rcorc is the cluster core radius (in pc), Vs (km s~^) 
the velocity scale, and (mo) the mean mass (in Mq) of all 
particles in thermal equilibrium in the central parts. 

Thus, significant mass segregation among the most mas- 
sive stars in the cluster core occurs on the local, central re- 
laxation time-scale (comparable to just a few crossing times, 
cf. Bonnell & Davies 1998), whereas a time-scale oc tr,h is 
required to affect a large fraction of the cluster mass. 

It should be kept in mind, however, that even the con- 
cept of a "local relaxation time" is only a general approxi- 
mation, as dynamical evolution is a continuing process. The 
time-scale for a cluster to lose all traces of its initial con- 
ditions also depends on the smoothness of its gravitational 
potential, i.e. the number of stars (Bonnell & Davies 1998: 
larger clusters are inherently smoother, and therefore mass 
segregation is slower than in smaller clusters with a grainier 
mass distribution), the degree of equipartition reached (e.g.. 
Hunter et al. 1995: full global, or even local, equipartition is 
never reached in a realistic star cluster, not even among the 
most massive species), and the slope of the MF (e.g., Light- 
man & Shapiro 1978, Inagaki & Saslaw 1985, Pryor, Smith 
& McClure 1986, Sosin 1997: flatter mass spectra will speed 
up the dynamical evolution, whereas steep mass spectra will 
tend to a higher degree of equipartition), among others. 

In addition, as the more massive stars move inwards 
towards the cluster centre, their dynamical evolution will 
speed up, and hence the dynamical relaxation time-scale for 
a specific massive species is hard to define properly. This 
process will be accelerated if there is no (full) equipartition 
(cf. Inagaki & Saslaw 1985), thus producing high-density 
cores very rapidly, where stellar encounters occur very fre- 
quently and binary formation is thought to be very effective 
(cf. Inagaki & Saslaw 1985, Elson et al. 1987b). In fact, the 
presence of binary stars may accelerate the mass segrega- 
tion significantly, since two-body encounters between bina- 
ries and between binaries and single stars are very efficient 
(e.g., Nemec & Harris 1987, De Marchi & Paresce 1996, Bon- 
nell & Davies 1998, Elson et al. 1998). This process will act 
on similar (or slightly shorter) time-scales as the conven- 
tional dynamical mass segregation (cf. Nemec & Harris 1987, 
Bonnell & Davies 1998, Elson et al. 1998). In summary, the 
time-scale for dynamical relaxation is a strong function of 
position within a cluster, and varies with its age. 



1.2 Primordial Mass Segregation 

Although a cluster will have lost all traces of its initial condi- 
tions on time-scales longer than its characteristic relaxation 
time, on shorter time-scales the observed stellar density dis- 
tribution is likely the result of dynamical relaxation and of 
the way that star formation has taken place. The process 
is in fact more complicated, as the high-mass stars evolve 
on the same time-scale as the lower mass stars (cf. Aarseth 
1999). In order to understand the process of mass segrega- 
tion in a cluster in detail, we have to get an idea of the 
amount of "primordial" mass segregation in the cluster. 

The nature and degree of primordial mass segregation 
is presumably determined by the properties of interactions 
of protostellar material during the star-forming episode in a 
cluster. In the classic picture of star formation (Shu, Adams 
& Lizano 1987), interactions are unimportant, and mass seg- 
regation does not occur. However, Fischer et al. (1998) con- 
clude that their observations of NGC 2157 seem to indicate 
the picture in which encounters at the early stages in a clus- 
ter's evolution enhance mass accretion due to the merging of 
protostellar clumps until the mass of these clumps exceeds 
the initial mass of a star to be formed. More massive stars are 
subject to more mergers, hence accrete even more mass (cf. 
Larson 1991, Bonnell et al. 2001a,b and references therein), 
and therefore dissipate more kinetic energy. In addition, they 
tend to form near the cluster centre, in the highest-density 
region, where the encounter-rate is highest (cf. Larson 1991, 
Bonnell et al. 1997, 1998, 2001a,b, BonneU & Davies 1998). 
This will lead to an observed position-dependent MF con- 
taining more low-mass stars at larger radii compared to the 
MF in the cluster centre (although low-mass stars are still 
present at small radii). This scenario is fully consistent with 
the idea that more massive stars tend to form in clumps and 
lower-mass stars form throughout the cluster (Hunter et al. 
1995, Brandl et al. 1996, and references therein). 

Although it has been claimed that the observed mass 
segregation in R136, the central cluster in the large star 
forming complex 30 Doradus in the LMC, is likely at least 
partially primordial (e.g., Malumuth & Heap 1994, Brandl 
et al. 1996) its age of ~ 3-4 Myr is sufficiently long for at 
least some dynamical mass segregation, in particular of the 
high- mass stars in the core (r < 0.5 pc), to have taken place 
(cf. Malumuth & Heap 1994, Hunter et al. 1995, Brandl et 
al. 1996). On the other hand, the presence of the high-mass 
Trapezium stars in the centre of the very young Orion Neb- 
ula Cluster (ONC; < 1 Myr, equivalent to ~ 3-5 crossing 
times; Bonnell & Davies 1998) is likely largely due to mass 
segregation at birth (Bonnell & Davies 1998, based on nu- 
merical simulations; Hillenbrand & Hartmann 1998, based 
on the appearance of the cluster as non-dynamically relaxed, 
and references therein). Bonnell & Davies (1998) show con- 
vincingly that the massive stars in the core of the ONC most 
likely originated within the inner 10-20% of the cluster. 

Hillenbrand & Hartmann (1998) argue that the young 
embedded clusters NGC 2024 and Monoceros R2 also show 
evidence for primordial mass segregation, since the outer 
regions of these clusters (and of the ONC as well) are not 
even one crossing time old. 
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2 THE DATA 

As part of Hubble Space Telescope (HST) programme GO- 
7307, we obtained deep WFPC2 V and /-band imaging of 
7 rich, compact star clusters in the LMC, covering a large 
age range. In de Grijs et al. (2001; Paper I) we presented 
the observational data for the two youngest clusters in our 
sample, NGC 1805 and NGC 1818, and discussed the depen- 
dence of the LPs on radius within each cluster. We found 
clear evidence for luminosity segregation within the inner 
~ 30" for both clusters, in the sense that the inner annular 
LPs showed a relative overabundance of bright stars with 
respect to the less luminous stellar population compared to 
the outer annular LPs. 

In this paper, we will extend our analysis to the associ- 
ated MPs and discuss the implications of our results in terms 
of the IMP and the star formation process. In Section ^ we 
will derive the MP slopes for both clusters, using a number 
of mass-to-luminosity (ML) conversions discussed in Section 
^ We will take care to only include main sequence stars be- 
longing to the clusters in our final MPs; to do so, we will ex- 
clude the field LMC red giant branch stars from the colour- 
magnitude diagrams (CMDs), with colours {V — I) > 0.67 if 
they are brighter than V — 22 mag (see the CMDs in John- 
son et al. 2001 for comparison). In fact, as can be seen from 
these CMDs, the true structure of the HR diagram above the 
main-sequence turn-off is very complex. Stellar populations 
of different masses overlap in colour-magnitude space, so 
that unambiguous mass determination from isochrone fits, 
for the handful to the few dozen stars populating these ar- 
eas in each cluster, is highly model-dependent (cf. Pig. 11 in 
Johnson et al. 2001). In a differential analysis such as pre- 
sented in this paper, the uncertainties involved in their mass 
determinations are too large and systematic (i.e., model de- 
pendent), so that we cannot include these stars in our anal- 
ysis. 

Table 1 of Paper I contains the fundamental parameters 
for our two young sample clusters. Por the analysis in this 
paper, however, we need to justify our choice for the adopted 
metallicity, age, and cluster mass in more detail. 

(i) Cluster metallicities - Por NGC 1805, metallicity 
determinations are scarce. Johnson et al. (2001) obtained an 
estimate of near solar metallicity, from fits to HST CMDs. 
The only other metallicity estimate available for NGC 1805, 
[Pe/H] ~ -0.30 (Mehani et al. 1994) is based on the average 
metallicity of the young LMC population and is therefore 
less certain. 

Abundance estimates for NGC 1818, on the other hand, 
are readily available, but exhibit a significant range. The 
most recent determination by Johnson et al. (2001), based 
on HST CMD fits, similarly suggests near-solar abundance, 
[Pe/H] « 0.0. Metallicity determinations based on stellar 
spectroscopy range from roughly [Pe/H] ~ —0.8 (Meliani et 
al. 1994, WiU et al. 1995, Oliva & Origha 1998) to [Pe/H] 
~ —0.4 (Jasniewicz & Thevenin 1994, Bonatto et al. 1995; 
see also Johnson et al. 2001). 

Por the purposes of the present paper, we will consider the 
cases of [Pe/H] = 0.0 and [Pe/H] = -0.5 for both clusters. 

(ii) Age estimates - Various age estimates exist for both 
clusters, which are all roughly consistent with each other, 
although based on independent diagnostics. The age range 
for NGC 1805 is approximately bracketed by logf(yr) = 



6.95-7.00 (cf. Bica et al. 1990, Barbaro & Olivi 1991, Santos 
Jr. et al. 1995, Cassatella et al. 1996) and logt(yr) = 7.6-7.7 
(cf. Barbaro & Olivi 1991), the most recent determinations 
favouring younger ages. We will therefore adopt an age for 
NGC 1805 of logi(yr) = T.Oto.l- 

Numerous age estimates are available for NGC 1818, on 
average indicating a slightly older age for this cluster than 
for NGC 1805. Most estimates bracket the age range be- 
tween t ~ 15 Myr (Bica et al. 1990, Bonatto et al. 1995, 
Santos Jr. et al. 1995, Cassatella et al. 1996) and t > 65 
Myr (cf. Barbaro & Olivi 1991), with the most recent esti- 
mates, based on HST CMD fits, favouring an age t ~ 20 — 30 
Myr (e.g., Cassatella et al. 1996, Grebel et al. 1997, Hunter 
et al. 1997, van Bever & Vanbeveren 1997, Pabregat & Tor- 
rejon 2000). We will therefore adopt an age of t ~ 25 Myr 
for NGC 1818, or logt(yr) = TAtH- 

(iii) Cluster masses - To obtain mass estimates for both 
clusters, we first obtained the total V-band luminosity for 
each cluster based on fits to the surface brightness profiles of 
our longest CPN exposures, in order to retain a sufficiently 
high signal-to-noise ratio even for the fainter underlying stel- 
lar component. We subsequently corrected these estimates 
of Lv,tot for the presence of a large number of saturated stars 
in the long CPN exposures (i.e, the observations where we 
located the cluster centre in the WFPC2/PC chip) by com- 
parison with the short CPN exposures. Although even in 
the short CPN exposures there are some saturated stars (cf. 
Section 3.3 in Paper I), their number is small (12 in NGC 
1805 and 18 in NGC 1818), so that we can get firm lower 
limits of logLv,tot{Lv,e) = 4.847 and 5.388 for NGC 1805 
and NGC 1818, respectively. 

Models of single-burst simple stellar populations (e.g., 
Bruzual & Chariot 1996), which are fairly good approxi- 
mations of coeval star clusters, predict mass-to-light (M/L) 
ratios as a function of age, which we can use to obtain 
photometric mass estimates for our clusters. This leads to 
mass estimates of Mtot = 2.8lo;8 x 10^ Mq (logMtot/M© = 
3.45to?5) for NGC 1805 and Mtot = 2.3lo J x 10'' Mq, or 
logMtot/Me = 4.35l°:J^ for NGC 1818. 

Our mass estimate for NGC 1805 is low compared to the 
only other available mass, Mtot = 6 x 10^ Mq (Johnson et al. 
2001), although their mass estimate, based on earlier simple 
stellar population models, is close to the upper mass allowed 
by our la uncertainty. 

Por NGC 1818, our mass estimate is entirely within the 
probable range derived by Plson, Pall & Preeman (1987), 
i.e. 4.1 < log Mtot /M0 < 5.7, depending on the M/L ratio, 
and Hunter et al.'s (1997) estimate of Mtot = 3 x W^Mq 
falls comfortably within our la uncertainty. Chrysovergis et 
al.'s (1989) determination of log Aftot/MQ = 4.69 is outside 
our la error bar; we speculate that the difference between 
our two estimates is due to a combination of the single- 
mass isotropic King cluster model used by them, versus our 
photometric mass determination, and a different treatment 
of the background stellar contribution. 



3 CONVERTING LUMINOSITY TO MASS 
FUNCTIONS 

The conversion of an observational LP (which we deter- 
mined for NGC 1805 and NGC 1818 in Paper I), in a given 
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passband i, (j){Mi), to its associated MF, ^(m), is not as 
straightforward as often assumed. The differential present- 
day stellar LF, dN/d<j){Mi), i.e. the number of stars in the 
absolute-magnitude interval [Mi, Mi + dMi], and the differ- 
ential present-day MF, dN/d^{m), i.e. the mass in the cor- 
responding mass interval [m,m + dm], are related through 
dN = -(t){M^)dMi = C(m)dm (Kroupa 2000), and therefore 

0(MO = -CM;^. (3) 

Thus, in order to convert an observational LF into a 
reliable MF, one needs to have an accurate knowledge of 
the appropriate ML - or mass-absolute-magnitude - rela- 
tion, dm/dMi. Empirical ML relations are hard to come by, 
and have so far only been obtained for solar-metallicity stars 
(e.g.. Popper 1980; Andersen 1991; Henry & McCarthy 1993, 
hereafter HM93; Kroupa, Tout & Gilmore 1993, hereafter 
KTG93). The ML relation is, however, a strong function of 
the stellar metallicity, and one needs to include corrections 
for hidden companion stars to avoid introducing a system- 
atic bias in the derived MF (e.g., KTG93, Kroupa 2000). 
For the conversion of the present-day MF to the IMF, one 
needs additional corrections for stellar evolution on and off 
the main sequence, including corrections for age, mass loss, 
spread in metallicity and evolution of rotational angular mo- 
mentum, or spin (cf. Scalo 1986, Kroupa 2000). Although 
the ML relation is relatively well-established for stars more 
massive than ~ O.8M0, our rather limited understanding 
of the lower-mass, more metal-poor stars, especially of the 
boundary conditions between the stellar interior and their 
atmospheres, have until recently severely limited the appli- 
cability of reliable ML relations to obtain robust MFs at the 
low-mass end. 

3.1 The mass-luminosity relation down to ~ 0.4AfQ 

As shown by Eq. (^), it is in fact the slope of the ML re- 
lation at a given absolute magnitude that determines the 
corresponding mass, which is therefore quite model depen- 
dent. This has been addressed in detail by, e.g., D'Antona 
& Mazzitelh (1983), Kroupa et al. (1990, 1993), Elson et al. 
(1995) and Kroupa & Tout (1997). 

The slope of the ML relation varies significantly with 
absolute magnitude, or mass. As shown by Kroupa et al. 
(1990, 1993) for solar-metallicity stars with masses m < 
IMq, it has a local maximum at My ~ 7, and reaches 
a minimum at My ~ 11.5 (see also Kroupa 2000). This 
pronounced minimum corresponds to a maximum in the 
present-day LF, while the local maximum at My ~ 7 cor- 
responds to the Wielen dip in the present-day LF of nearby 
stars (e.g., Kroupa et al. 1990, D'Antona & Mazzitelli 1996, 
and references therein). 

The local maximum in the derivative of the ML relation 
at Mv ~ 7 (m ~ 0.7^/©) is caused by the increased impor- 
tance of the E[~ opacity in low-mass stars with decreasing 
mass (KTG93, Kroupa & Tout 1997). 

The ML relation steepens near AIv ~ 10 (m ~ 0.4 — 
O.5M0), due to the increased importance of H[2 formation in 
the outer shells of main sequence stars, which in turn leads 
to core contraction (e.g., Chabrier & Baraffe 1997, Baraffe 
et al. 1998, Kroupa 2000). 

Given the non-linear shape of the ML relation and the 



small slope at the low-mass end, any attempt to model the 
ML relation by either a polynomial fit or a power-law de- 
pendence will yield intrinsically unreliable MFs (cf. Elson et 
al. 1995, Chabrier & Mera 1997), in particular in the low- 
mass regime. This model dependence is clearly illustrated 
by, e.g., Ferraro et al. (1997), who compared the MFs for 
the GC NGC 6752 derived from a variety of different ML 
relations at that time available in the literature. 

3.2 Age and metallicity dependence and 
corrections for binarity 

The exact shape of the ML relation is sensitive to metallic- 
ity; metallicity changes affect the stellar spectral energy dis- 
tribution and therefore the (absolute) magnitude in a given 
optical passband (cf. Brewer et al. 1993). In fact, it has been 
argued (cf. Baraffe et al. 1998) that, although the V-band 
ML relation is strongly metallicity-dependent, the K-hand 
ML relation is only a very weak function of metal abun- 
dance, yielding similar _R'-band fluxes for [AI/H] = —0.5 
and [Af/jy] = 0.0. Although the ML relation is currently 
relatively well-determined for solar-metallicity stars with 
m > O.8Af0, at low metallicities the relation remains very 
uncertain. This is partially due to the lack of an empirical 
comparison, and to our still relatively poor understanding 
of the physical properties of these stars, although major ef- 
forts are currently under way to alleviate this latter problem 
(e.g., the recent work by the Lyon group). 

Fortunately, as long as we only consider unevolved main 
sequence stars, age effects are negligible and can therefore 
be ignored (cf. Brewer et al. 1993, Ferraro et al. 1997). This 
applies to the current study for the stellar mass range con- 
sidered. 

Finally, stellar populations contain in general at least 
50 per cent of multiple systems. The immediate effect of 
neglecting a significant fraction of binary stars in our LF- 
to-MF conversion will be an underestimate of the resulting 
MF slope (cf. Fischer et al. 1998, Kroupa 2000). We wiU 
come back to this point in Section ^.3| . 

3.3 A comparison of mass-luminosity relations 

3. 3. 1 Luminosity-to-mass conversion for stars with 
masses below IMq 

Since no empirical ML relations are available for low-mass, 
low-metallicity main sequence stars, a test of the goodness 
of ML relations in this regime must therefore bear on the 
comparison of different models. Several recent studies have 
adopted this approach (e.g., Alexander et al. 1997, Ferraro 
et al. 1997, Kroupa & Tout 1997, Piotto, Cool & King 1997, 
Saviane et al. 1998). 

For solar-metallicity stars in the mass range 0.1 < m < 
IMq, Leggett et al. (1996) and Kroupa & Tout (1997) con- 
cluded that, although all models considered provided rea- 
sonable fits to the empirical ML relation, the Baraffe et al. 
(1995) theoretical ML relations provided the best overall 
agreement with all recent observational constraints. On the 
other hand, Bedin et al. (2001) show that these are poor at 
low metallicity. It should be noted that the Baraffe et al. 
(1995) models were based on grey model atmospheres. 

Both Piotto et al. (1997) and Saviane et al. (1998), from 
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Figure 1. Empirical anti theoretical ML relations, (a) - solid bullets: HM93; open circles: Andersen (1991). (b) - solid line: HM93 fit; 
dotted line: KTG93 and Kroupa & Tout (1997) semi-empirical ML relation; dashed line: Chabrier et al. (1996) theoretical ML relation 
for m < 0.6Mq, based on a third-order polynomial fit. (c) ~ Theoretical ML relations for subsolar abundances: Alexander et al. (1997; 
solid lines, for [M/H] = —1.3, —1.5, —2.0 [top to bottom]), and Baraffe et al. (1997; dashed line, [M/H] = —1.5). For comparison, the 
solar-abundance ML relation of Baraffe et al. (1997) is also shown (dash-dotted line), (d) - Observational data for m > I.OMq stars 
(Andersen 1991), and - for m < 5Mq - theoretical models by GBBCOO for solar abundance (solid line) and [M/H] = —1.3 (dashed line). 



a comparison of largely the same theoretical ML relations 
available in the literature with observational data for the 
low-metallicity Galactic GCs NGC 6397 ([Fe/H] ~ -1.9) 
and NGC 1851 ([Fe/H] ~ -L3), respectively, concluded that 
the Alexander et al. (1997) theoretical ML relations for the 
appropriate metallicity provided the best match for masses 
J; 0.6 — O.8A'/0. Similar conclusions were drawn by Piotto 



et al. (1997) for three other Galactic GCs, M15, M30 and 
M92. Alexander et al. (1997) themselves found a good to ex- 
cellent overall agreement between their models and those of 
the Lyon group, in particular the updated Chabrier, Baraffe 
& Plez (1996) ones, which employ the most recent non-grey 
model atmospheres. 

Figure ^ shows the available empirical data, on which 
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these comparisons are based for solar-metallicity stellar pop- 
ulations. The filled bullets represent the HM93 sample; the 
open circles the higher-mass Andersen (1991) binary stars. 
In panel (b), we show the m < 2Mq subsample. Overplotted 
are the best-fitting relation of HM93 (solid line), the fit to 
their semi-empirical ML relation (dotted line) of KTG93 and 
Kroupa & Tout (1997), and the theoretical ML relation of 
Chabrier et al. (1996; dashed line) for 0.075 <m< 0.6Mq. 
The figure shows that the observational data allow for signif- 
icant local differences in the slope of the solar-metallicity ML 
relation; these uncertainties propagate through the deriva- 
tive of the relation when converting LFs to MFs. 

The theoretical ML relation for solar abundance by 
Chabrier et al. (1996) closely follows the most recent semi- 
empirical ML relation compiled by Kroupa (KTG93, Kroupa 
&Tout 1997). In Fig.|l|c, we compare the current theoretical 
ML relations for subsolar metallicity: the solid lines repre- 
sent the Alexander et al. (1997) ML relations for (top to bot- 
tom) [M/H] ~ —1.3, —1.5, and —2.0^; for reasons of clarity, 
we only show the [M/H] = —1.5 ML relation of Baraffe et 
al. (1997), but the spread due to metallicity differences is 
similar to that shown by the Alexander et al. (1997) rela- 
tions. The most significant differences between both sets of 
models are seen at masses m > OAMq. This is likely due 
to the slightly different treatment of the stellar atmospheres 
and radiative opacities. Finally, for comparison we also show 
the solar-metallicity theoretical ML relation of Chabrier et 
al. (1996) and Baraffe et al. (1997). 

We plot the derivatives of the ML relations as a func- 
tion of absolute visual magnitude in Fig. ^ From Fig. ^ it is 
immediately clear that the empirical fit to the HM93 ML re- 
lation inherently leads to unreliable luminosity-to-mass con- 
versions because of the two sharp discontinuities in the slope. 

Figure ^ shows the metallicity dependence of the slope 
of the ML relation; the solid lines represent the Alexander 
et al. (1997) ML relations with [M/H] = -1.3,-1.5, and 
-2.0, peaking from right to left. The Baraffe et al. (1997) 
models (cf. the dashed line, for [M/H] = —1.5) closely follow 
the Alexander et al. (1997) ones. For comparison, we have 
also included the solar-abundance model of Chabrier et al. 
(1996) and Baraffe et al. (1997), as in panel (a). 



3.3.2 The more massive stellar population 

The main uncertainties for the luminosity evolution of stars 
with masses m > O.SMq are in the treatment of the degree 
of mass loss and convective core overshooting. Girardi et al. 
(2000, hereafter GBBCOO) and Girardi (2001, priv. comm.) 
computed a grid of stellar evolutionary models for stars in 

t We applied the procedure outlined in Ryan & Norris (1991) 
to convert the Alexander et al. (1997) metallicities (Z) to [M/H] 
values: [M/H] ^ [O/H] = [O/Fe] -|- [Fe/H], with [O/Fe] = -1-0.35 
for [Fe/H] < -1, and [O/Fe] = -0.35x [Fe/H] for -1 < [Fe/H] 
< 0. Therefore, a stellar population with an observed [Fe/H] = 
-1.5 corresponds to a model with [Z] = log( Z / Zq ) = [M/H] ~ 
— 1.15, and Z ~ 1.35 X 10""^. This procedure has been shown to 
apply to halo subdwarfs; these oxygen-enriched abundances are 
characteristic of old stellar populations in the Milky Way (cf. 
Baraffe et al. 1995). Gilmore & Wyse (1991) have shown that the 
element ratios in the LMC are significantly different, however. 



the mass range 0.15 < m < 7Mq for metallicities between 
^ and 1.5 times solar, using updated input physics, as well 
as moderate core overshooting. 

In Fig. |l|d, we show the observational ML relation of 
Andersen (1991) for stars with masses m > l.OM©. In addi- 
tion, we have plotted GBBCOO's (2000) theoretical ML rela- 
tions for solar metallicity (solid line) and for [M/H] = — 1.3 
(dashed line) for stars less massive than ~ 5Mq. These mod- 
els include moderate core overshooting, but the presence or 
absence of this process in the models does not significantly 
change the resulting ML relation for this mass range. 

3.3.3 Comparison for our young LMC clusters 

Based on the comparison and discussion in the previous sec- 
tions, for the conversion of our observational (individual) 
stellar magnitudes (Paper I) to masses, and thence to MFs 
we will use 

• the empirical ML relation of IIM93. This ML relation 
is defined for stars with My > 1.45. 

• the KTG93 and Kroupa & Tout (1997) semi-empirical 
ML relation for stars with Mv > 2.00, with an extension to 
AIv ~ —3 by adoption of Scalo's (1986) mass-Mv relation. 

• the parametrisation of these by Tout et al. (1996, 
hereafter TPEH96), valid for masses in the range —1 < 
log m/M© < 2; we converted the corresponding bolometric 
luminosities to absolute IZ-band magnitudes using the bolo- 
metric corrections of Lejeune, Cuisinier & Buser (1998). The 
TPEH96 ML relations are given as a function of metallicity 
from Z = 10"* to Z = 0.03. 

• the GBBCOO models. For solar metallicities, the models 
for their youngest isochrone of 60 Myr are defined for stars 
with —3.381 < Mv < 12.911, while for the subsolar abun- 
dance of Z = 0.008 this corresponds to -4.832 < Mv < 
12.562. 

Although we argued that the Baraffe et al. (1998, 
hereafter BCAH98) models employ the most recent input 
physics, their mass range, m < I.OMq, precludes us from 
using their models, since completeness generally drops be- 
low our 50% limit for m < O.8i\f0, thus leaving us with too 
few data points for a useful comparison. 

In Fig. ^ we compare the mass estimates based 
on HM93, KTG93, BCAH98 and GBBCOO to TPEH96's 
parametrisation, for both solar and subsolar ([Fe/H] — —0.5) 
metallicities and ages of ~ 10 and 25 Myr. Significant differ- 
ences are appreciated among the individual models, in par- 
ticular between TPEH96 on the one hand and the high-mass 
end (log tu/Mq > 0.3) of KTG93 and between TPEH96 and 
the models of BCAH98, which show systematic deviations 
from the one-to-one relation indicated by the dashed line. 
It is therefore not unlikely that the differences among the 
models dominate the uncertainties in the derived MF and - 
ultimately - in the IMF slope (see also Bedin et al. 2001, who 
reached a similar conclusion in their analysis of the Galactic 
GC M4). We will quantify these effects in the next section. 

In Fig. ^ we show examples of the derived global cluster 
MFs, normalized to unit area, for three of the ML conver- 
sions adopted in this paper, KTG93, TPEH96 and GBBCOO. 
All MFs are shown for solar metallicity and for mass bins 
corresponding to luminosity ranges that exceed our 50% 
completeness limits (see Paper I). The MFs, and all other 
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Figure 2. The slope of the ML relation, (a) - Solar-metallieity ML relations: HM93 (dashed line, empirical), KTG93 and Kroupa & 
Tout (1997) (dotted line, semi-empirical), and Chabrier et al. (1996) and BarafTe et al. (1997) (solid line, theoretical), (b) - Theoretical 
ML relations for subsolar abundances: Alexander et al. (1997; solid lines, for [M/H] = —1.3, —1.5, —2.0 [peaking from right to left]), and 
Baraffe et al. (1997; dashed line, [M/H] = —1.5). For comparison, the solar-abundance ML relation of Chabrier et al. (1996) and Baraffe 
et al. (1997) is also shown (dash-dotted line, as in (a)). 



MFs discussed in this paper, were corrected for incomplete- 
ness and background contamination following identical pro- 
cedures as for the LFs in Paper I. Significant systematic 
effects result from the adoption of any given ML conver- 
sion, as can be clearly seen. For comparison, we also plot 
a fiducial Salpeter IMF, which appears to be a reasonable 
approximation for the global cluster MFs in the mass range 
(-0.15 < logm/Mo < 0.8). 



4 QUANTIFICATION OF MASS 
SEGREGATION EFFECTS 

4.1 Radial dependence of luminosity and mass 
functions 

The effects of mass segregation can be quantified using 
a variety of methods. The most popular and straightfor- 
ward diagnostic for mass segregation effects is undoubtedly 
the dependence of MF slope, F — A log ^(m)/ Am (where 
f (m) oc mF), on cluster radius. 

In Fig. 1^ we plot the derived MF slopes as a function 
of cluster radius for our four adopted ML conversions and 
assuming three different fitting ranges in mass: 

• —0.15 < log m/M0 < 0.30 for all conversions; 

• -0.15 < logm/Mo < 0.70 for KTG93, TPEH96 and 
GBBCOO; and 

• -0.15 < logm/A/o < 0.85 for TPEH96 and GBBCOO. 

The adopted radial ranges for our annular MFs are indi- 
cated by the horizontal bars at the bottom of each panel. Al- 



though we see a large spread among models and mass fitting 
ranges, clear mass segregation is observed in both clusters 
at radii r < 20", well outside the cluster core radii (indi- 
cated by the vertical dotted lines). We have also indicated 
the Salpeter (1955) IMF slope, F = -1.35 (dashed horizon- 
tal lines) . While both the global MFs (which are dominated 
by the inner, mass-segregated stellar population) and the 
annular MFs near the core radii appear to be consistent 
with the Salpeter IMF slope (cf. Vesperini & Heggie 1997; 
for the cluster models assumed here, Rcovc ~ Rh), the MFs 
beyond the cluster radii where mass segregation is signifi- 
cant (i.e., r > 30")are characterised by steeper slopes, i.e., 
relatively more low-mass stars compared to high mass stars 
than found in the inner cluster regions. This result holds for 
all mass fitting ranges, all ML conversions considered, solar 
and subsolar ([Fe/H] — —0.5) metallicity, and both clusters. 

The error bars in Fig. ^ represent the formal uncertainty 
in the fits; the systematic uncertainties are clearly greater, 
and a strong function both of the adopted ML conversion 
and of the mass range used for the fitting of the MF slopes. 
In all cases, the larger mass ranges used for the fitting result 
in steeper MF slopes than the smaller (lowest-mass) range, 
thus presenting clear evidence for non-power law shaped 
MFs. In addition, the TPEH96 parametrisation results in 
systematically steeper slopes for all cluster-metallicity com- 
binations. 

Because of the strong model dependence, in particu- 
lar because of the sensitivity to the choice of ML relation, 
and the accuracy of the corrections for incompleteness and 
background star contamination of single power law fits to 
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logf mass / Mq) [Tout et al. 1996] 



Figure 3. Comparison of tlie mass estimates for NGC 1805 (open circles) and NGC 1818 (filled circles) resulting from various ML 
relations. We used TPEH96's parametrisation as comparison ML relation because of its large mass range. 



© 2001 RAS, MNRAS 000, 000-000 



Mass segregation in young LMC clusters II. 9 




Figure 4. Examples of the derived global cluster MPs, normalized to unit area, for three of the ML conversions adopted in this paper, 
KTG93 (solid lines), TPEH96 (dashed lines) and GBBCOO (dotted lines). All MPs are shown for solar metallicity and for mass bins 
corresponding to luminosity ranges that exceed our 50% completeness limits (see text). The dash-dotted lines show the Salpeter IMF 
corresponding to the mass range under consideration. 



the annular MFs, in Fig. ^ we introduce a more robust 
characterisation of the presence of mass segregation in these 
two young star clusters. We quantified the deviations of the 
high-mass range of the annular MFs from the global MF 
following a similar procedure as defined in Paper I for the 
LFs. All annular MFs were normalized to the global MF 
in the range 0.00 < log m/M© < 0.20, where the effects of 
mass segregation - if any - are negligible, as shown above; 
subsequently, we determined the sum of the differences be- 
tween the global and the scaled annular MFs in the common 
mass range 0.20 < logm/M© < 0.60, E(Ar). We adopted 
log m/M© = 0.60 as upper mass limit, so that we could 
compare the results of all three ML conversions. 

From Fig. ^ it follows that both clusters are mass seg- 
regated within R ~ 30"; beyond, the deviations become rel- 
atively constant with increasing radius. It is also clear (i) 
that there is no appreciable difference between the strength 
of the mass segregation in the two clusters, and (ii) that the 
scatter among the data points from the different models is 
relatively small. We therefore conclude that we have been 
able to quantify the effects of mass segregation in a fairly 
robust way by thus minimising the effects of the choice of 
ML conversion. 

It is most likely that the effect referred to as mass segre- 
gation is indeed due to a positional dependence of the ratio 
of high-mass to low-mass stars within the clusters, and not 
to different age distributions ("age segregation"). Johnson et 
al. (2001) have shown that for the high-mass stars in both 
clusters the LF is indeed just a smoothed (almost) coeval 
CMD, for which age effects are only of second order im- 



portance. It is possible that for very low stellar masses, i.e. 
pre-main-sequence stars, age segregation may play a more 
important role, but this applies only to stars well below the 
50% completeness limits for both clusters. 

Finally, we point out that it is generally preferred to use 
star counts rather than surface brightness profiles to mea- 
sure mass segregation effects (e.g., Elson et al. 1987b, Cher- 
noff & Weinberg 1990). Elson et al. (1987b) argue that the 
use of surface brightness profiles by themselves, although ini- 
tially used to study mass segregation (e.g., da Costa 1982, 
Richer & Fahlman 1989), is limited in the sense that one 
cannot distinguish between these effects and any significant 
degree of radial velocity anisotropy in a cluster's outer re- 
gions (see also Chernoff & Weinberg 1990). In addition, a 
comparison between results obtained from star counts and 
from surface brightness profiles does not necessarily trace 
the same system, since star counts are generally dominated 
by main-sequence (and subgiant-branch) stars, while surface 
brightness profiles mostly trace the giants (and also subgiant 
stars) in a cluster (cf. Elson et al. 1987b). 



4.2 Core radii 

Conclusive results on the presence of mass segregation in 
clusters can also be obtained by examining the core radii 
of specific massive stellar species (e.g., Brandl et al. 1996), 
or - inversely - by measuring the mean stellar mass within 
a given radius (e.g., Bonnell & Davies 1998, Hillenbrand & 
Hartmann 1998). However, it may not always be feasible to 
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Figure 5. Mass function slopes as a function of cluster radius. The symbols are coded as follows: stars — HM93, squares — KTG93, 
triangles — TPEH96, circles — GBBCOO; open symbols correspond to a mass range —0.15 < logm/MQ < 0.30, filled symbols to 
—0.15 < log m/M0 < 0.70, and crossed open symbols extend up to logm/MQ = 0.85. The vertical dotted lines indicate the cluster core 
radii; the Salpeter slope is represented by the horizontal dashed lines. The radial ranges over which the MF slopes were determined are 
shown by horizontal bars at the bottom of each panel. The right-hand subpanels show the overall MF slopes for the clusters as a whole; 
the data points are spread out radially for display purposes. 



use this diagnostic, since the individual stellar masses of the 
cluster members need to be known accurately, thus provid- 
ing an additional observational challenge. In addition, the 
results depend critically on which stars are used to obtain 
the mean mass, and can be severely aflected by small num- 
ber statistics (cf. Bonnell & Davies 1998). 

Figure ^ and Table ^ show the dependence of the de- 



rived cluster core radius on the adopted magnitude (or mass) 
range. Core radii were derived based on fits to stellar num- 
ber counts - corrected for the effects of incompleteness|^ and 

i We only used magnitude (mass) ranges for which the complete- 
ness fractions, as determined in Paper I, were at least 50%. 
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Figure 6. Deviations of the annular MFs from the global MF as a function of radius, as discussed in the text. Filled symbols were 
obtained using the ML relation of KTG93, open symbols are based on TPEH96's ML conversion, and crossed open symbols result from 
the GBBCOO models. The horizontal bars at the bottom of the figure indicate the radial range used to obtain the data points, the cluster 
core radii are indicated by the vertical dotted lines. 
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Figure 7. Core radii as a function of magnitude (mass). The filled circles are the core radii after correction for the effects of 
(in)coinplctencss, area covered by the observations, and background stars; the open circles are not background subtracted and serve 
to indicate the uncertainties due to background correction. We have also indicated the mean cluster core radii, obtained from surface 
brightness profile fits (dotted lines). The horizontal bars at the bottom of the panels indicate the magnitude ranges used to obtain the 
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Table 1. Cluster core radii as a function of mass 



Magnitude 
range {V) 


log m, [ A/q 
(central) 


NGC 
(") 


1805 
(pc) 


NGC 
(") 


1818 
(pc) 


15.0 - 17.0 


0.90 


7.30 


1.85 


11.65 


2.95 


17.0 - 18.0 


0.75 


5.74 


1.45 


17.45 


4.41 


18.0 - 19.0 


0.60 


5.34 


1.35 


13.61 


3.44 


19.0 - 20.0 


0.44 


5.27 


1.33 


15.01 


3.80 


20.0 - 21.0 


0.30 


13.97 


3.53 


18.07 


4.57 


21.0 - 22.0 


0.16 


19.39 


4.91 


30.74 


7.78 


22.0 - 23.0 


0.09 


21.09 


5.34 


31.63 


8.00 


23.0 - 24.0 


0.00 


23.86 


6.04 


29.55 


7.48 


15.0 - 24.0 




16.96 


4.29 


23.92 


6.05 



Note: The range (15.0 ^ V < 24.0) represents the cluster average, corrected 
for completeness, area covered by the observations, and the background 
population. 

background contamination (cf. Paper I) - of the generalised 
fitting function proposed by Elson et al. (1987a), in the lin- 
ear regime: 

MO=Mo||l+(^0'j , (4) 

where ^(r) and /lo are the radial and central surface bright- 
ness, respectively, 7 corresponds to the profile slope in the 
outer regions of the cluster, and J?corc ~ a(2^/^ — 1)^/^ « iih- 
Equation reduces to a modified Hubble law for 7 = 2, 
which is a good approximation to the canonical King model 
for GCs (King 1966). 

For both NGC 1805 and NGC 1818 we clearly see the ef- 
fects of mass segregation for stars with masses logm/MQ > 
0.2 {Mv < 2.4; m > I.6M0). It is also clear that the bright- 
est four magnitude ranges, i.e. masses log m/M© > 0.4(m > 
2.5M0), show a similar concentration, while a trend of in- 
creasing core radius with decreasing mass (increasing mag- 
nitude) is apparent for lower masses. The larger scatter for 
NGC 1818 is due to the smaller number of stars in each 
magnitude bin compared to NGC 1805; for NGC 1818 the 
associated uncertainties are determined by a combination of 
the scatter in the derived core radii and background effects, 
while the uncertainties for NGC 1805 are dominated by the 
effects of background subtraction. 

We have also indicated the core radii obtained from 
profile fits to the overall surface brightness profiles of the 
clusters. It is clear that these are dominated by the mass- 
segregated high-mass (bright) stars. 



5 A COMPARISON OF MASS FUNCTION 
SLOPES - TRACING THE IMF? 

5.1 Comparison with previously published results 

Few studies have published MPs of sufficient detail and 
quality for the two young LMC clusters analysed in this 
paper to allow useful comparisons. Santiago et al. (2001) 
published global MFs and MFs determined in the annulus 
4.9 < 7? < 7.3 pc (19.4 < R < 28.8") for both NGC 1805 
and NGC 1818, based on the same observations used for 
the present study, while Hunter et al. (1997) published the 
global MF as well as the core MF derived from the HST 
WFPC2/PC chip of NGC 1818. In Fig. | we compare our 



results with those of Hunter et al. (1997) and Santiago et al. 
(2001). 

Although Santiago et al. (2001) used a different ML 
conversion than done in the present paper, the slopes and 
the dependence on the adopted mass fitting range they de- 
rived for their annular MFs are fully consistent with the 
range seen at this radius, using any of the ML relations 
employed by us. However, their and Hunter et al.'s (1997) 
global MF slopes are somewhat shallower than ours. The 
difference is sufficiently small, however, that it can be ex- 
plained as due to the combination of different ML relations 
and a different treatment of the background stellar popula- 
tion (see Paper I for a discussion of the latter). Hunter et 
al. (1997) found no significant difference in MF slope be- 
tween the core MF (F = -1.21 ± 0.10) and the global MF 
(r — —1.25 ± 0.08), in the mass range between 0.85 and 
9Mo(-0.07 < logm/A'/0 < 0.95). The PC field of view 
samples the inner r ~ 18"; we have also included their core 
data point in Fig. |[ Although Hunter et al.'s (1997) PC 
MF slope appears to be slightly steeper than most of our 
MF slope determinations at these radii, this can easily be 
explained as due to a combination of the uncertainties in the 
ML conversion used (as evidenced by the range in MF slopes 
seen in Fig. ^) and the intrinsic curvature of the MF. From 
Fig. ^it follows that the MF slopes get steeper if increasingly 
higher mass stars are included in the fitting range (compare 
the location of the open, filled and open-crossed symbols, 
which indicate increasing mass fitting ranges). Hunter et 
al.'s (1997) PC MF slope determination is based on a mass 
range extending up to 9Af0, while our largest fitting range 
only includes stars < 7Mq. The observed curvature in the 
MF will therefore result in a slightly steeper slope for the 
Hunter et al. (1997) MF slope, although still fairly similar to 
the slopes determined using our greatest mass fitting ranges 
at these radii (cf. Fig. |^. 

With regard to the use of different ML relations, it is 
worth noting here that the difference between the MF slopes 
derived by us and those of both Hunter et al. (1997) and 
Santiago et al. (2001) may be largely due to the adopted 
isochrones: while for older clusters the MF slopes for main- 
sequence stars are almost independent of age, small dif- 
ferences between MF slopes as a function of age are ap- 
preciated for younger stellar populations. This difference is 
in the sense that using isochrones for older stellar popula- 
tions will result in slightly shallower MF slopes. This is the 
most likely explanation for the slight shift between our MF 
slopes (based on solar neighbourhood-type stellar popula- 
tions) and those of Hunter et al. (1997) and Santiago et al. 
(2001), who both used younger isochrones to obtain their 
mass estimates. However, the expected steepening in MF 
slope from evolved to young stellar populations is almost 
entirely contained within the observed spread in MF slope 
in both clusters (cf. Hunter et al. 1997). In fact, for the 
ML conversion based on the GBBCOO models, we used their 
youngest isochrone (at t = 6 x 10^ yr); the resulting MFs 
are shown as circles in Fig. ^. Hunter et al. (1997) showed 
that if they had used a 30-40 Myr isochrone instead of the 
one at 20 Myr used by these authors, would make the IMF 
slope appear shallower hy AT ~ 0.15. Thus, it appears that 
the GBBCOO slopes are entirely consistent with the slopes 
obtained from the other ML relations, however, which em- 
phasizes our statement that any slope difference due to the 
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Figure 8. Comparison of MF slopes. The small symbols represent our NGC 1805 and NGC 1818 data points; the symbol and line coding 
is as in Fig. ^ The large circles represent Santiago et al.'s (2001) MF slopes; the large squares those of Hunter et al. (1997) for NGC 
1818, while the symbol coding is as for the small symbols. The radial ranges are again indicated by the horizontal bars in each panel. 
The arrow with vertical error bar in the panel of NGC 1818 (located at r = 70") shows the global MF slope derived by Hunter et al. 
(1997) and its associated la uncertainty. Solar abundances were assumed. 



adoption of different isochrones is accounted for by the ob- 
served spread in MF slopes. 

A further comparison is provided by the resuh of Will 
et al. (1995), who obtained a MF slope of F = -1.1 ± 0.3 
for stars between 2 and 8 Mq (or logm/M© = 0.30 - 0.90). 
This slope is significantly shallower than any of the slopes 
derived by us, in particular in view of the observed steepen- 
ing of the MF slope when including increasingly higher-mass 
stars (Section |4.1[ ). A similar discrepancy was already noted 
by Hunter et al. (1997), who argued that this difference (and 
their greater uncertainties) was most likely due to the diffi- 
culty of obtaining reliable photometry from crowded ground- 
based images, while they were not able to resolve stars in the 
cluster centre, nor detect stars of similarly faint magnitudes 
(low masses) as possible with HST observations. 



5.2 Primordial or dynamical mass segregation? 

The key question is whether the observed mass segregation 
in both young LMC star clusters is the result of the process 
of star formation itself or due to dynamical relaxation. In 
Fig. ^ we have plotted the half-mass relaxation time as a 
function of mass, using Eq. (|l]) and the mass-dependent core 
radii of Fig. ^ For comparison, we have also indicated the 
ages of both clusters. For NGC 1805 significant dynamical 
mass segregation is expected to have occurred out to its 
half-mass radius for stars more massive than about 3 Mq 
(logm/M© ~ 0.48), while for NGC 1818 this corresponds to 
stars exceeding ~ 8Mq {\ogm/ Mq > 0.90). However, from 
Fig. ^ it follows that mass segregation becomes significant 
for masses m > 2.5Mq, out to at least 20 — 30", or 3 — 6i?corc. 




• NGC 1B05 
O NGC 1818 



2 4 6 8 

Mass (Mq) 

Figure 9. Half-mass relaxation time as a function of mass for 
NGC 1805 and NGC 1818. The best age estimates for both clus- 
ters are indicated by horizontal dashed lines. 

Dynamical mass segregation in the cluster cores will 
have occurred on 10 — 20 x shorter time-scales, in particular 
for the more massive stars (cf. Eq. (^)). In fact, if the cluster 
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contains a significant amount of gas, e.g., Mgas > Mstars (cf. 
Lada 1991, Bonnell & Davies 1998), this will increase the 
cluster's gravitational potential, and thus the virialised stel- 
lar velocity dispersion (Bonnell & Davies 1998). Therefore, 
in this case a larger number of two-body encounters, and 
hence time, is required to reach a dynamically relaxed state. 
Thus, the relaxation time estimates obtained by consider- 
ing only the contributions of the cluster's stellar component 
should be considered lower limits, especially for young star 
clusters, which are generally rich in gas. 

Elson et al. (1987b) estimated the central velocity dis- 
persion in NGC 1818 to be in the range 1.1 < ctq < 6.8 
km s~^. Combining this central velocity dispersion, the core 
radius of ~ 2.6 pc, and the cluster age of ~ 25 Myr, we esti- 
mate that the cluster core is between ~ 5 and ~ 30 crossing 
times old, so that dynamical mass segregation in the core 
should be well under way. Although we do not have velocity 
dispersion information for NGC 1805, it is particularly in- 
teresting to extend this analysis to this younger ('^ 10 Myr) 
cluster. We know that its core radius is roughly half that of 
NGC 1818, and its mass is a factor of ^ 10 smaller. Simple 
scaling of Eq. (|l]) shows then that the half-mass relaxation 
time of NGC 1805 is ~ 4-5x as short as that of NGC 1818; 
if we substitute the scaling laws into Eq. (^), we estimate 
that the central velocity dispersion in NGC 1805 is > lOx 
smaller than that in NGC 1818. From this argument it fol- 
lows that the cluster core of NGC 1805 is < 3 — 4 crossing 
times old. 

However, since strong mass segregation is observed out 
to ~ 6i?core and ~ 37?core iu NGC 1805 and NGC 1818, 
respectively, for stellar masses in excess of ~ 2.5M0, it is 
most likely that significant primordial mass segregation was 
present in both clusters, particularly in NGC 1805. Although 
this was initially suggested by Santiago et al. (2001), we have 
now substantiated this claim quantitatively. Relevant to this 
discussion is the study by Bonnell & Davies (1998), who 
found that whenever a system of massive stars is found at 
the centre of a young star cluster, like the Trapezium stars in 
the ONC, a major fraction of it most likely originated in the 
inner parts of the cluster. N-body simulations are currently 
being carried out to investigate the fraction of massive stars, 
and their mass range, that will have to have originated in 
the cluster centres to result in the observed distribution. We 
will include these in a subsequent paper (de Grijs et al., in 
prep.; Paper III). 

5.3 The slope of the cluster mass function 
We will now return to the discussion of Figs. |5| and ^. In 



section 



4.1 



we showed (i) that the slope of the global cluster 
MF is relatively well approximated by that at the cluster 
core radius; (ii) that at the cluster core radius the effects 
of strong mass segregation are still clearly visible; and (Hi) 
that in the outer cluster regions, the slope of the (annular) 
cluster MFs approaches a constant value. 

Within the uncertainties, we cannot claim that the 
slopes of the outer MFs in NGC 1805 and NGC 1818 are 
significantly different. Starting with the work by Pryor et 
al. (1986) and McClure et al. (1986), it is expected that 
clusters with similar metallicities exhibit similar MF slopes. 
However, Santiago et al. (2001) claimed to have detected a 
significantly different MF slope for both the global and their 



annular MF between both clusters. As is clear from Fig. g|, 
their annular MFs were not taken at sufficiently large radii 
to avoid the effects of mass segregation. Since both clusters 
are affected by mass segregation in a slightly different way 
(which may just be a reflection of the difference in their 
dynamical ages), it is not surprising that Santiago et al.'s 
(2001) annular MFs exhibit different slopes. 

Recent studies show that the actual value of the MF 
slope may vary substantially from one region to another, 
depending on parameters such as the recent star formation 
rate, metallicity, and mass range (cf. Brandl et al. 1996). 
The outer cluster regions of R136/30Dor (Malumuth & Heap 
1994, Brandl et al. 1996), M5 (Richer & Fahlman 1987), M15 
(Sosin & King 1997) and M30 (Sosin 1997) are all charac- 
terized by MF slopes F < 2.0. BonneU et al. (2001b) ex- 
plain this rather steep MF slope naturally as due to the 
process of star formation and accretion itself, resulting from 
the combination of a gas dominated and a stellar dominated 
regime within the forming cluster. This results in a double 
power law IMF, where the lower mass stars have a shallower 
slope and the high-mass IMF slope is steeper (F « —1.5 and 
~2 < r < —2.5, respectively), due to the different accretion 
physics operating in each regime (i.e., tidal-lobe accretion 
versus Bondi-Hoyle accretion for low-mass and high-mass 
stars, resp.). 

if there is a significant per- 



As mentioned in section 3.2 



centage of binary or multiple stars in a star cluster, this 
will lead us to underestimate the MF slope. In other words, 
the MF slopes determined in this paper are lower limits, 
because we have assumed that all stars detected in both 
clusters are single stars. Elson et al. (1998) found a fraction 
of ~ 35 ± 5 per cent of roughly similar mass binaries (with 
mass of the primary ~ 2 — 5.5M0) in the centre of NGC 
1818, decreasing to ~ 20 ± 5 per cent in the outer regions of 
the cluster, which they showed to be consistent with dynam- 
ical mass segregation. It is not straightforward to correct the 
observed LFs for the presence of binaries, in particular since 
the binary fraction as a function of brightness is difficult to 
determine. If we consider the information at hand, (i) the 
small (~ 15%) gradient in the total binary fraction derived 
by Elson et al. (1998) for the inner ~ 3 core radii of NGC 
1818; (a) the result of Rubenstein & Bailyn (1999) that the 
binary fraction increases at fainter magnitudes in the Galac- 
tic GC NGC 65752; and (in) the similarity of our annular 
LFs for faint magnitudes for all annuli and both clusters, we 
conclude that the effects of a binary population in either or 
both of NGC 1805 and NGC 1818 are likely smaUer than 
the observed differences in LF shapes as a function of radial 
distance from the cluster centres. This is corroborated by 
Elson et al.'s (1998) result for NGC 1818. 



6 SUMMARY AND CONCLUSIONS 

We have reviewed the complications involved in the conver- 
sion of observational LFs into robust MFs, which we have 
illustrated using a a number of recently published ML re- 
lations. These ML relations were subsequently applied to 
convert the observed LFs of NGC 1805 and NGC 1818, the 
two youngest star clusters in our HST programme of rich 
compact LMC star clusters, into MFs. 

The radial dependence of the MF slopes indicate clear 



© 2001 RAS, MNRAS 000, 000-000 



Mass segregation in young LMC clusters II. 15 



mass segregation in both clusters at radii r < 20 — 30", well 
outside the cluster core radii. This result does not depend 
on the mass range used to fit the slopes or the metallicity 
assumed. In all cases, the larger mass ranges used for the 
fitting result in steeper MF slopes than the smallest mass 
range dominated by the lowest-mass stars, thus presenting 
clear evidence for non-power law shaped MFs. Within the 
uncertainties, we cannot claim that the slopes of the outer 
MFs in NGC 1805 and NGC 1818 are significantly different. 
Wc also argue that our results are consistent with previously 
published results for these clusters if we properly take the 
large uncertainties in the conversion of LFs to MFs into 
account. The MF slopes obtained in this paper are in fact 
lower limits if there is a significant fraction of binary stars 
present in the clusters. 

The global cluster MFs (which axe dominated by the 
inner, mass-segregated stellar population) and the annular 
MFs near the core radii appear to be characterised by sim- 
ilar slopes, the MFs beyond the cluster radii where mass 
segregation is significant (i.e., r > 30") are characterised by 
steeper slopes. It is, however, not unusual for star clusters 
to be characterised by rather steep MF slopes; Bonnell et al. 
(2001b) explain this naturally as due to the different accre- 
tion physics operating in the low and high-mass star forming 
regime. 

We analysed the dependence of the cluster core radius 
on the adopted magnitude (mass) range. For both clusters 
we clearly detect the effects of mass segregation for stars 
with masses logm/M0 > 0.2 (m > I.QMq). It is also clear 
that stars with masses log m/MQ > 0.4(m > 2.5M0) show a 
similar concentration, while a trend of increasing core radius 
with decreasing mass (increasing magnitude) is apparent for 
lower masses. The characteristic cluster core radii, obtained 
from profile fits to the overall surface brightness profiles, are 
dominated by the mass-segregated high-mass stars. 

We estimate that the NGC 1818 cluster core is between 
~ 5 and ~ 30 crossing times old, so that dynamical mass seg- 
regation in its core should be well under way. Although we 
do not have velocity dispersion information for NGC 1805, 
by applying scaling laws we conclude that its core is likely 
< 3 — 4 crossing times old. However, since strong mass seg- 
regation is observed out to ~ 6i?coro and ~ 3-Rcorc in NGC 
1805 and NGC 1818, respectively, for stellar masses in ex- 
cess of ^ 2.5Mq, it is most likely that significant primordial 
mass segregation was present in both clusters, particularly 
in NGC 1805. We are currently investigating this further 
using N-body simulations. 
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